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A master equation with a Lindblad structure is derived, which describes the interaction of a 
test particle with a macroscopic system and is expressed in terms of the operator valued dynamic 
structure factor of the system. In the case of a free Fermi or Bose gas the result is evaluated in the 
Brownian limit, thus obtaining a single generator master equation for the description of quantum 
Brownian motion in which the correction due to quantum statistics is explicitly calculated. The 
friction coefficients for Boltzmann and Bose or Fermi statistics are compared. 
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The study of the dynamics of a particle coupled to a general many-body system plays a relevant role in modern 
quantum physics, both with respect to foundations and applications of quantum theory. On the one hand it provides 
a most simple example of quantum dynamics of a non isolated system, possibly offering a manageable arena for a 
truly microscopic approach, which might shed some light on mechanisms of dissipation and decoherence [Q: these 
\ last two issues are now of outstanding relevance in connection with the rapidly growing experimental ability to deal 
ON ■ with thoroughly quantum mechanical phenomena, checking for their coherence properties (at the single particle level 
think for example of the recent cavity QED and ion trapping experiments , while at many-body level Bose- Einstein 
condensation is a most interesting example Q). On the other hand plenty of interesting physical problems may be 
modeled in this way and among these, in particular, motion or diffusion of charged or neutral particles in gases, liquids 
or solids. The interaction of a test particle with a dilute or noninteracting gas is strictly connected to the problem of 
a quantum generalization of the Boltzmann equation, whose everlasting relevance has recently been stressed by the 
experimental realization of quantum degenerate samples of weakly interacting bosons or fermions j^Q] : in fact for the 
study of these systems resort is often made to a quantum Boltzmann transport equation J5| . A particularly interesting 
situation arises if the mass M of the test particle is much bigger than the mass m of the particles which make up 
the gas: the so called Brownian motion, which serves as a paradigmatic example in the description of irreversible and 
dissipative processes. The description of the phenomenon is still debated at a quantum level (see || and references 
cited therein), even though well settled by now at classical level in terms of Langevin or Fokker-Planck equations 
(it took however almost a century from the observation by Brown to the first successful theoretical description by 
Einstein, which led to the first example of fluctuation-dissipation relation, linking friction and diffusion coefficients). 
A class of models, usually named quantum Brownian motion Q, given by time evolutions with a non Hamiltonian part 
mapping the algebra of operators at most bilinear in the operators x and p of the particle into itself, seem to be the 
most natural candidate in order to obtain equations of motion analogous to the classical ones, leading in particular to 
a friction force proportional to velocity. On the mathematical side generators of time evolution semigroups satisfying 
these requirements have been fully characterized through the property of complete positivity, which formally amounts 
to the requirement that positivity of the time evolution is preserved even in presence of coupling without interaction 
to another system and leads to a typical expression for the generators of these semigroups, known as Lindblad 
structure || . This has led to a wide literature developing this axiomatic approach || , together with a large number of 
more or less phenomenological models in which similar structures are obtained, though not always preserving complete 
positivity (in this connection see fi"c|| ). Though warranting positivity of the statistical operator, complete positivity 
is by itself no fundamental requirement as recently stressed [ [Tl"| , so that despite its extensive use in many fields of 
physics, ranging from quantum optics to quantum communication, the study of the conditions and approximations 
under which it emerges from microphysical models is strongly desirable. 

In recent work the derivation at a fundamental level of a completely positive master equation for a Brownian 
particle interacting with a free Boltzmann gas has been given, based on a microphysical model developed for the 
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description of particle matter interaction |12j,|6[|. The Lindblad equation thus obtained can be written with a single 
generator and temperature dependent friction and diffusion coefficients were determined in terms of the scattering 
cross section. In this paper we give a major extension of the previous model, keeping also quantum statistics of the 
gas into account. Moreover, before going over to the Brownian limit, in which the ratio between the masses is much 
smaller than one, one sees that the generator of the master equation is expressed in terms of the dynamic structure 
factor of the medium, first introduced by van Hove jl3f| - This turns out to be true also for an interacting system, thus 
linking in full generality the dynamics of the test particle to the density fluctuations of the system [see Eq. (Jigj)] . The 
property of complete positivity is retained in the general case under some requirements on the energy dependence 
of the dynamic structure factor, which are exactly fulfilled in the case of Boltzmann particles dealt with in || . The 
Brownian limit is then considered, thus obtaining the correction at finite temperature due to quantum statistics to 
the equation describing quantum Brownian motion [see Eq. (|2q)]. In terms of the fugacity z this correction takes a 
remarkably simple form [see Eq. ( 

The paper is organized as follows: in Sec. II we consider the general structure of the master equation and its 
connection to the dynamic structure factor; in Sec. Ill we obtain the correction due to quantum statistics to the 
master equation describing quantum Brownian motion, together with the relationship between the friction coefficients 
for Boltzmann or quantum statistics; in Sec. IV we comment on our results indicating potential future developments. 
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II. GENERAL STRUCTURE OF THE MASTER EQUATION IN TERMS OF THE DYNAMIC 

STRUCTURE FACTOR 

Let us briefly recall the structure of the master equation obtained in [Q for the description of the subdynamics of 
a particle interacting with a macroscopic system supposed to be at equilibrium. The result is valid on a time scale 
r much longer than microphysical collision time and describes an interaction through two-particle collisions given by 
the full T matrix. The master equation is given by: 

f = 4 [Ho, 8] + \Y. l^U HL]*L*, §}] , (1) 

with 
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where Ho is the Hamiltonian for the particle and g its statistical operator, while g m = X^^ICXCI i s the statistical 
operator for matter at equilibrium, ir^ being the statistical weights related to its spectral decomposition. The vectors 
| A) and |£) are eigenvectors of the macrosystem Hamiltonian H m with eigenvalues E\ and E^ respectively, similarly 
|fc) and \h) denote eigenvectors of Ho with eigenvalues Ek, and Eh- In writing the equation we have neglected the 
slow energy dependence of the T matrix, which would have brought a commutator term proportional to the forward 
scattering amplitude, diagonal in momentum representation. The terms other than the commutator in are linked 
to the dissipative behavior, which cannot be obtained in a Hamiltonian formalism. Interactions at microphysical level 
are typically translationally invariant, so that a general Ansatz for the T matrix describing two-body interactions is 
given by T\ = J d 3 x J d 3 y (x)u* k (y)t(x — y)uh(y)ip(x), where ip 1 , ip are field operators for the macrosystem. We 
now consider a homogeneous system, so as to use as quantum numbers momentum eigenvalues, thus obtaining with 
a Fourier transform an expression depending only on the modulus of the momentum transfer: 



T h = 5 Pr,+ Pk -P h + P A\Pn -P v \) b v b u> ( 2 ) 
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where b\ b denote creation and destruction operators in the Fock space of the macrosystem. Restricting to the 
case of a free gas of Bose or Fermi particles, the eigenvectors of H m can be characterized as a set of occupation 
numbers n a relative to particles with a given momentum p a , so that |£) = |{rtf }), and the matrix element (A|6* fr^l^) 
can be readily evaluated restricted to the primed sum for A ^ £, since in the case A = £ the contributions to the 
master equation (|l|) cancel out. Denoting by q = p^ — p n the momentum transferred to the test particle and by 

AE fiq (p) — + ^ 2 J* 2M ~ 1m tne difference in energy before and after the collision (M being the mass of 

the test particle with momentum p, m the mass of the gas particles), and supposing the statistical operator g to be 
quasi-diagonal in momentum representation, according to its slow variability, one sees that (Q) for a free test particle 
reduces to || 
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accordingly, z denoting the fugacity, determined by the requirement J2fi( n n) = N, and (3 — l/(kBT) the inverse 
temperature. It is worthwhile introducing the more compact notation 
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where n denotes the density of particles in the gas and the function S B/F is in fact positive definite. The integral in 
(0) can be explicitly calculated both for bosons and fermions giving at finite temperature the result 
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where <r(q,p) = ^ [(1 + a)q 2 + 2a(p ■ q)] is expressed in terms of the dimensionless variable a = m/M, giving the 
ratio between the masses. Expression (H) is exactly the dynamic structure factor for a free Bose or Fermi gas at finite 
temperature, as one could also directly realize from (|i|) ]l4| or from the equivalent expression in terms of momentum 
transfer q and energy transfer E = + ^ (note that we use as variables momentum and energy transferred to the 
particle) 
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where the dependence on the transferred momentum is in this case actually only through the modulus. In the following 
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we will use, according to convenience, both notations S(q,p) and S(q, E), where E = AE q (p) — E p+q —E p — 2M + ^- 
is the energy transfer. The dynamic structure factor is an important physical quantity of direct experimental access, 
essentially depending on the statistical properties of the macrosystem and the kinematics of the collision, appearing 
in the expression of the inelastic differential cross section for a particle interacting with a macroscopic sample. The 
relation between differential cross section and dynamic structure factor was first derived by van Hove in the case of 
neutron scattering Jl3| and for scattering from state p to state p' = p + q is given by 
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The dynamic structure factor is expressed in the general case as Fourier transform of the time dependent pair 
correlation function with respect to energy and momentum transfer, according to 



S(q,E) = J^i JdtJ d 3 xe^ Et -^ J d 3 y (N(y)N(x + y,t)) , 



(7) 



where N(y) denotes the local particle density for the macroscopic system and (. . .) the ensemble average. Alternatively 
the dynamic structure factor may be written in terms of the Fourier transform of the density operator N(y), given by 
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thus obtaining 
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Expression (Q) through the relation (g) allows a determination of the equilibrium fluctuations of the system in terms 
of scattering experiments fl4| (think for example of the very interesting applications in the case of neutron scattering 
from different states and isotopes of Helium |l5|]). Coming back to (|5|) we note that in the limit of very small fugacity 
z <C 1 one recovers the result for Maxwell Boltzmann particles 
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which in terms of momentum and energy transfer may also be written 
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Recalling expression (Q) for the dynamic structure factor one immediately realizes that the master equation given 
in (^) can be written in terms of the dynamic structure factor and exactly exhibits a Lindblad structure provided 
the dynamic structure factor evaluated at the arithmetic mean of p and p' equals the geometric mean of its values 
at the two points. This identity holds true without approximations in the case of expression (j^]) for a Boltzmann 
gas in the Brownian limit considered in 0. In the general case this factorization relies on an approximation linked 
to the quasi-diagonality of the statistical operator. Keeping the linear relation between E and p into account, the 
approximation necessary in order to retain complete positivity can be most meaningfully written 
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and will depend on the smoothness of the energy dependence of S in the relevant energy region (note that the neglected 
terms are at least quadratic in the energy difference). Exploiting ([tl]) Eq. (||) can be cast in the following Lindblad 
structure granting positivity of the time evolution 
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which may be also written in a more manifest Lindblad form 
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introducing the following generator depending on the operators x and p 
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This is a remarkably simple result since L B/F only depends on the generator of translations in momentum space and 
the operator valued dynamic structure factor. Let us note that equation ( |l2| ) or equivalently ( p"3| ) is invariant under 

translation and rotation and in particular a statistical operator of the canonical form g oc e^^^w is a stationary 
solution. If instead of a free gas one considers a more general medium characterized by a dynamic structure factor 
S(q,p), provided the interaction between particle and medium still satisfies translation invariance as in (0) and an 
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approximation of the form (|Tl|) holds, the master equation (||) still has the form ( |l2| ) or equivalently ( |l3| ) with L(q, p, x) 
given by 



L(q,p,x) = e^y/S(q^j 

and therefore retains a completely positive structure. To prove this we go back to which in the case of a 
homogeneous system using (^J) can be written 

f = -i[H ,£] + ^) (15) 
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We now introduce the momentum transfer q = — p rl , q' = p , — p , and the Fourier transform p q of the density 
operator given by (|J), so that relabelling the indexes (|l^) becomes 
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where the primed sum over q and q' means that the contribution for q = q' = is left out, since in this case the two 
terms cancel out. To proceed further we express the denominators in terms of a Laplace transform, according to 
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where (...) denotes the ensemble average over g m , p q (t) the Heisenberg operator e + i H ™ lt p q e~% H ' nt , and the more 
compact notation AE q (p) — E p+q — E p for the energy transfer has been used. Since the system is supposed to be 
homogeneous, the correlation function selects the contributions for which q = q' , and exploiting the identity 
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We can now meaningfully undo the Laplace transform, coming to 
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If we now exploit the quasi-diagonality of linked to its slow variability, thus substituting in the denominators of 
the first term p, p' with the symmetric expression i(p + p'), we obtain the expression 
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The correlation functions appearing in ( |18[ ) are exactly the dynamic structure factor multiplied by N and evaluated 
for a momentum transfer q and energy transfers AE q ( p+ 2 p ) and AE q (p) respectively, as can be seen by comparison 
with (Q). To see under which conditions the obtained master equation ( |l5| ) takes a Lindblad structure we consider 
an approximation of the form (^) , which will generally depend on the smoothness of the energy dependence of the 
dynamic structure factor, but is actually less demanding than it might seem, since in the expression ( |l8| ) one has to 
consider a sum over p and p' with the matrix elements (p\g\p') of the statistical operator. In the continuum limit we 
therefore obtain the master equation 
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which still has the form fll2|), but is much more general since now the dynamic structure factor does not necessarily 
describe a free gas. This result allows for the extension of the usefulness of the master equation to cases in which the 
correlation function cannot be directly evaluated, but a suitable phenomenological model is available, e.g., determined 
in terms of scattering experiments. 



III. QUANTUM BROWNIAN MOTION AND QUANTUM STATISTICS 

We are now interested in the Brownian limit a — m/M <C 1, considering the dynamics of a free particle interacting 
through collisions with a gas of much lighter particles. Having an expression valid for both a Fermi or Bose gas it 
is particularly interesting to evaluate the correction brought about by quantum statistics to the typical models of 
quantum Brownian motion. In the limit a <C 1 expressions (^) and ([To]) become respectively 
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still satisfying the principle of detailed balance In the Boltzmann case, as mentioned above expression (|2l| ) 

exactly fulfills ([n]) and the generator in jl^ ) takes the particularly simple form L B/F (q, p,x) cx es'^e^w?'^ so that 
one obtains for an isotropic medium the master equation given in H 
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To recover the equation describing quantum Brownian motion one goes over to small momentum transfer, strongly 
favored by the kinematics of the collisions, considering terms up to second order in q or equivalently bilinear in x 
and p, thus obtaining an equation in close analogy to the classical description, with a friction force proportional to 
velocity. The result for a Boltzmann gas is 
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and has the particular feature that it can be written in Lindblad form in terms of a single generator Q . Starting from 
( pp| ) one can perform the same limit of small momentum transfer corresponding through the physical interpretation 
of the dynamic structure factor to the macroscopic, long wavelength properties of the system, thus calculating the 
correction due to quantum statistics to the master equation describing quantum Brownian motion. To do this one 
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considers the Taylor expansion of the logarithms in (g0|) , leading to the following compact expression as a power series 
in the fugacity z 
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which has to be substituted in (|l2|), keeping terms up to second order in q. The result one obtains is actually 
remarkably simple: the operator structure is not changed, nor the simple generator feature, but the fugacity appears 
through the expression z/(l =F z) rather than linearly. For a Bose or Fermi gas at finite temperature one has 
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where the coefficient z/(l =p z) at finite temperature is actually well defined because z is in the range < z < 1 for 
Bose particles and positive for Fermi particles. 

Eq. ( p6[ ) expressing the correction due to quantum statistics in the equation describing quantum Brownian motion, 
together with Eq. (TL9) giving a completely positive time evolution for a particle interacting with a macroscopic system 
at equilibrium in terms of a momentum displacement operator and the dynamic structure factor of the system, are 
the main results of this paper. Comparing ( p6j ) with ( p^ ) one sees that the friction coefficient given in the Boltzmann 
case by 
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enhanced or suppressed according to statistics. Both ( p3| ) and (|2^) retain the property of complete positivity satisfied 

by (|l3|), are invariant under translation and rotation and admit a stationary solution of the canonical form g oc e - ^^ . 
The single generator feature is due to the fact that the coefficients satisfy the relationship D pp D xx = fi 2r j 2 /4. 



IV. SUMMARY AND OUTLOOK 



We have considered the problem of the motion of a test particle interacting through collisions with a fluid, following 
the approach outlined in [|l2|,[l6],[r7j , which has already been successfully applied to the case of neutron optics 1 1^] . The 
microscopic derivation allows some insights into the conditions under which a master equation of the Lindb ad type, 
driving a completely positive time evolution, can be obtained, thus giving a concrete physical example contributing 
to the debate on the relevance of complete positivity jn]]. Provided the statistical operator is sufficiently diagonal 
in momentum representation with respect to the energy dependence of the dynamic structure factor, the master 
equation ( |l9| ) is obtained, where only quantities of physical interest appear: the scattering cross section for the 
single two-body collisions, given by the square of the T matrix; the generator of translations in momentum space 
and the dynamic structure factor, keeping the statistical properties of the medium into account, combined through 
the expression L(q, p,x) = e^ q ' x y / S(q, p). This structure is remarkably simple and describes a dynamics in which 
the motion of the test particle is linked through this particular two-point correlation function to the spectrum of 
spontaneous fluctuations of the system. Starting from this general structure and explicitly calculating the dynamic 
structure factor for the case of a free gas one can consider the particularly relevant case of Brownian motion, when the 
test particle is much heavier than the particles making up the gas. In the case of a Boltzmann gas one recovers, for 
small momentum transfer, a typical structure of generator of quantum Brownian motion, given by Eq. (p3[), in which 
all coefficients are determined and the dissipative part of the generator depends linearly on the fugacity. The case 
of a quantum gas is also considered, and in this case the generator has the structure (|26j), with the dissipative part 
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depending on the fugacity z through the expression z/(l =p z), thus giving the connection ( p8[ ) between the friction 
coefBcients in the different cases. 

We hope that this fundamental study on the general features of a master equation describing the motion of a test 
particle in a gas, putting in major evidence the dynamic structure factor and showing the relationship between this 
structure and the equation, analogous to the Fokker- Planck equation, describing quantum Brownian motion, could 
be a sound starting point for future extensions and applications, especially in connection with degenerate regimes at 
very low temperatures, where the dynamic structure factor is now being intensively studied both at experimental jnj 
and theoretical level potl. 
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